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ABSTRACT 

A fine analysis of the photometric standard star 
29 Psc has been performed using flux constant hydrogen- 
line-blanketed model atmospheres. 

The observable quantities, Balmer and Paschen slopes, 
Balmer discontinuity and H and Hel lines were used to 
determine the final model parameters: T ^ = 15150°K, 

log g =4.08 and N(He)/N(H) = 0.027. 

Lines from 14 elements besides H and He were 
identified, five elements (C, Mg, Si, Ti and Fe) had a 
sufficient number of unblended lines with known gf values 
to attempt an abundance determination. 
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Visiting Astronomer, Kitt Peak National Observatory, 
which is operated by the Association of Universities 
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National Science Foundation. 


I . INTRODUCTION 


The star 29 Psc (HD 224926) was chosen by Oke (1964) 
as a photoelectric standard because it was within 17° of the 
celestial equator, relatively bright (m = 5.1), earlier 
than A1 and not known to be variable or abnormal. Bless 
(1969) pointed out that 29 Psc had a small Balmer disconti- 
nuity for its spectral type (B8 III) from which it was 
concluded that the star was similar to helium-rich stars. 
However, in order to fit all the observable quantities it 
will be shown that the star is really helium-weak and it 
is hotter than B8, approximately B6 . This paper describes 
the atmospheric and abundance analyses of 29 Psc. The 
models used were hydrogen-line-blanketed models (Klinglesmith 
1971) . 

Section II describes the observations that were used; 
Section III describes the theoretical calculations used to 
compare with the observations and Section IV describes the 
fitting procedure used to determine the final model parameters. 
Sections V and VI describe the abundance and mass-radius 
determination. The rotational velocity determination and 
profile fits are discussed in Section VII. 

II. THE OBSERVATIONS 
a) Spectra 

Four coude-spectra were obtained at the Kitt Peak 
National Observatory 84-inch telescope on January 25 and 26, 
1970. The spectra were taken on IlaO plates and developed 
with the necessary calibration plates in D-76 for 12 minutes. 
Table 1 lists the pertinent data for the exposures. The 
linear dispersion was 8.9 S/mm with a resolving power of 
approximately 25,000. The useful spectral range was \3700 
to \4800. The density was recorded on magnetic tape with a 


Joyce-Lobel microdensitometer and converted to intensity in 
the computer. The four separate spectra were averaged to 
form one composite spectrum. It was possible to measure 
equivalent widths as small as 40 milliangstroms . Table 2 
lists the observed lines with their identifications and 
equivalent widths. Lines for which there are no equivalent 
width measurements were considered weak enough to be seen 
but not measured. 

b) Continuous Energy Distribution 

Since 29 Psc is a photoelectric standard, there exists 
good continuum data in the literature (Oke 1964 and Hayes 
1970) . The two sets of data are essentially in agreement 
(Hayes, Oke and Schild 1970). 

Figure 1 shows the original data from Hayes (1970) and 
the renormalization for the current calibration of Vega (Oke 
and Schild 1970) . It can be seen that the revision has 
the largest effect on the Balmer slope, in fact, changing 
the sign of the Balmer slope from negative to positive. This 
permits a consistent temperature to be determined from the 
3 continuum observables as will be described in Section IVa. 

c) Colors 

Crawford ( 1963) obtained UBV and HB colors for 29 Psc: 

U-B = -0.50, B-V = -0.3, HB = 2.708, E (U _ B ) =0.0 and =0.0. 
From its position in the color-color and the He vs. (U-B) q 
plots (Figure 1 and 2 of Crawford’s paper) it can be inferred 
that 20 Psc is indeed a B6 not a B8 star. Irwin (1961) used 
29 Psc as a comparison standard for his Southern Cepheid 
Photometry and did not find any indication of variation in 
the colors. 
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III. THEORETICAL CALCULATIONS 


The model stellar atmospheres used in the analysis were 
flux constant, hydrogen line-blanketed models (Klinglesmith 
1971). The usual assumptions of hydrostatic, radiative and 
thermodynamic equilibrium were made. The sources of opacity 
included were the bound-bound transitions of the Lyman and 
Balmer series of hydrogen, the bound-free and the free-free 
transitions of hydrogen, helium and their ions both positive 
and negative. 

The line profiles and equivalent widths of H H. and 

Y i o 

H were computed using the "ESW" formulation of Stark broadening 
and the Hel lines were computed using Stark broadening data 
from 2 different sources: for \4121 and 4713 from Griem (1964) 

and for \4026, 4388 and 4471 from Barnard, Cooper and Shamey 
(1969) . 

The line absorption coefficient for the metal lines 
used in the abundance analysis was a Voight function. The 
damping constants used are listed in Table 3 along with the 
other atomic data . In each case it was attempted to use the 
most realistic value that was available. 

The theoretical Balmer discontinuity was defined as 


D 3647 * - 2 5 lo *10 (> v ( 3647 + >/F v ( 3647-)) 


where log F^ (3647 + ) was obtained by linear extrapolation using 
log (4546) and log F (4217) . The rather peculiar wavelengths 
result from the quadrature points used in the model atmosphere 

calculations; \4217 is halfway between H. and H and X4586 is 

6 Y 

halfway between and HB. 

The Balmer and Paschen slopes are defined as 



-2.5 log 10 (F v (3647")F v (3421 + ))/AX(p) ; A\ 


0.0226 p 


- -2.5 log 1Q GM 5508) /F v ( 5050))/ A\( u ); AX = 0.0458 u 

Again, the peculiar wavelengths are a result of the quadrature 
points used in the model calculations. 


4 


IV. FITTING THE OBSERVABLES 


The observable quantities: Balmer and Paschen slopes, 

Balmer discontinuity, and equivalent widths of the hydrogen 
and helium lines were measured in a manner identical to the 
way they were computed from the theoretical model atmospheres. 
The theoretical quantities were plotted as a function of 
effective temperature for each log g with one curve for each 
helium to hydrogen ratio. A plot of the T 0ff - log g pairs 
that fit each of the observables represents the locus of 
the models that would give the observed value; again, one 
curve for each abundance. The loci have different slopes 
in the T 0ff - log g plane due to the different T 0ff - log g 
dependence of each observable. However, the loci of a given 
observable for different helium to hydrogen ratios are 
essentially parallel. The intersection of two or more loci 
for a given helium to hydrogen ratio determines a model that 
fits those two or more observables. Figure 2 shows these 
loci for 29 Psc. The individual loci will be described in 
the following sections. 

a) The Continuum Observables 

The values for the continuum observables for the revised 
photoelectric energy distribution of 29 Psc (Figure 1) are 

°3647 = -0 • 71 » = 1 - 80 and ^p " 2 - 68 - 

Each one of the 3 continuum observables gave a similar 
loci, i.e. no gravity dependence and very little temperature 
dependence for a given helium to hydrogen ratio, however, there 
is a temperature spread of about 1000°K when the helium to 
hydrogen ratio varies from 0 to 1.5. Hence, the loci of each 
observable at one helium to hydrogen ratio were averaged 
together. Figure 2 shows the resultant "averaged loci".. They 
are labeled H-continuum. 


b) Hydrogen Lines 


The equivalent widths of the hydrogen lines on the 
other hand do show a dependence on both T gff and log g as 
well as abundance. Three lines were used (Hy, H g and H e ) 
and their individual loci were averaged for each helium to 
hydrogen ratio. At a given log g the total spread in 
temperature between the loci of H y , H g and H e was 800°K 
or less. The average of the 3 loci were plotted in Figure 
2 - labeled H-lines. 

At this point it is possible to draw a line (marked HLC 
in Figure 2) that shows the interceptions of the hydrogen 
line loci with the hydrogen continuum loci. This line 
indicates the parameters of the models that would best fit 
the hydrogen observables. Notice that the helium to hydrogen 
rttio would vary along the line. 

c) The Helium Lines 

The equivalent widths of 5 Hel lines for which good 
Stark broadening data exist were also used to obtain "averaged 
loci" at four helium to hydrogen ratios (0.013, 0.028, 0.125 
and 1.5). The equivalent widths for the two smallest ratios 
were obtained by using a pure hydrogen atmosphere and assuming 
the helium was an impurity. The five lines were 4026, 4388, 
4471, 4121 and 4713. The individual loci at a given helium 
to hydrogen ratio had a spread of ~ 1500°K at a given log g. 

The dashed line marked HHE is the line of intercepts 
between the hydrogen lines and the helium lines. The inter- 
section of the lines HHE and HLC at T gff = 15150°K and log g - 
4.08 determine the effective temperature and log g for 29 Psc. 
In order to determine the helium to hydrogen ratio at this 
T ef f and log g, three models were computed with the previously 
determined T gff and log g and different helium to hydrogen 
ratios (0.013, 0.028 and 0.044). The profiles and equivalent 
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widths of the 5 Hel lines were recomputed for each model. 

The variation of with N(He)/N(H) is shown in Figure 3. 

The crosses indicate the observed for each of the lines. 

The individual N(He)/N(H) and these averages are given in 
Table 3. Thus, from the hydrogen continuum, the hydrogen 
lines and the helium lines a final model with the parameters 
T eff - 15150 + 700; log g - 4.08 + 0.15 and N(He/N(H) - 0.027 
± 0.012 has been determined. 

V. ABUNDANCE ANALYSIS 

Using the same techniques that have been used in previous 
analyses (Hunger and Klinglesmith 1969 and Klinglesmith , 

Hunger f Bless and Millis 1970) the abundances of five elements 
were determined. The only new feature of analysis is the 
inclusion of variable partition functions (Sparks and Fischel 
197] ) . Table 4 lists the lines used in the abundance 
analysis, their atomic data and the derived mass fractions. 
Table 5 gives the relative number fractions for 29 Psc and 
a normalized log N i for 29 Psc, E stars (Aller 1961) and 
the cosmic abundances (Allen 1963). From this, we can see 
that only the carbon abundance appears normal. Magnesium 
is underabundent (only one line was used) whereas Si, Ti and 
F e are overabundant by about an order of magnitude or more. 

The spread in the individual determination for Ti and Fe 
lines indicates that the gf values and/or the damping constants 
are still not well enough known. 

VI. MASS AND RADIUS 

There is no direct distance measure for 29 Psc, however, 
two indirect methods can be used in order to derive an estimate 
of the mass and radius. The spectral type that corresponds 
to 15000°K is B6 and the luminosity that corresponds to log g - 
4 is V or IV (Blaauw 1963). This implies a mean M y of -1.1 
and since there is no absorption (Crawford 1963) a distance of 
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« 175 parsecs results. Using Petrie's (1958a, b) H y vs. 
luminosity relation one obtains M y « -0.8 and a distance of 
150 parsecs. The bolometric correction for the final model 
is -1.32 which gives an average M Bq1 - -2.32. The mass and 
radius can be determined from the relations (see Klemola, 1961) 
5 log R - 42.4-10 log T eff - M B(j1 

and 

log M = log g + 2 log R -4.38. 

This results in M pa 7.4 M and R pa 3.8 R . These values 

© © 

are consistent with the mass and radius of a normal B star 
and should not be interpreted any more strongly than that. 

VII. ROTATIONAL VELOCITY 

In order to attempt profile fits it is necessary to know 
the rotational velocity or rather v r sin i. For 29 Psc this 
was determined from the Mgll 4481 line. First, the abundance 
of Mg was determined from the equivalent width. Then, with 
that abundance, an unrotated profile was computed and rotated 
using the expression given by Pecker and Schatzman (1959) , 


+1 

P r ' I_ L p u 0 { l 


v„ sin i 
r 

c 



tt/2(1-x 2 ) 
1 + 


i 1 2 

+ 2(1 -x^) 

-2 

3 p 


dx. 


This expression describes the effect on the unrotated profile, 

P , on a spherical star with a linear limb-darkening law, where 
0 is the limb-darkening coefficient. 

The Mgll 4481 line was rotated at 50, 60 and 70 km/sec 
and the best fit was obtained with 60 km/sec. 

Using these same values the Hel and hydrogen lines were 
also rotated, Figure 4 shows the results. For the Hel lines 
3 values of N(He)/N(H) are shown at each v r sin i. The best 


fit is obtained for a KHe)/N(H) « 0.028 and V r sin i - 60 km/sec. 
For the H lines only the 60 km/sec profiles are shown. The 
two calculated profiles represent N(He)/N(H) of 0.013 and 
0.125. The fit is good in the wings but the observed cores 
are still too deep. 

The metal lines whose equivalent widths were greater 
than 0.l8 were also rotated and again 60 km/sec was the 
v r sin i that was obtained if the unrotated profile had the 
abundance determined from the equivalent width measurement. 

VIII. CONCLUSIONS 

The standard LTE type of model atmosphere abundance 
analysis has shown that the observed data for 29 Psc can 
be consistently fitted by a model in which helium is under- 
abundant by a factor of 4 compared to the current estimate 
for N(He)/N(H) = 0.106 (Leckrone 1971) for normal main-sequence 
stars. Also the overabundance of the F e group elements along 
with the underabundance of helium are consistent with the 
idea that 29 Psc is a hot Ap star. 
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FIGURE CAPTIONS 


Fig. 1 - The emergent flux for 29 Psc . The dote are Hayes 
original data, the crosses are Hayes data corrected for 
tlva new calibration of Vega jy Oke and Schild. 

Fig. 2 - Atmospheric parameter determination; see text for 
description. The final model has T e ^ • 15150 and log 
g = 4.08. 

Fig. 3 - The equivalent widths of the Hel lines as a function 
of the Helium to Hydrogen ratio. The crosses are the 
observed values of W for 29 Psc. 

A 

Fig. 4 - Helium and Hydrogen line profile fits for 29 Psc. 

For the 5 Hel lines, 3 values of N(He)/N(H) are plotted 
0.013( ), 0.028( ) and 0.043(**O. The computed pro- 

files have been rotated to the value shown over the 
curves. The broad line is the observed curve. The 
hydrogen lines are rotated to 60 km/sec and shown for 2 
N(Hc) /N(H) (0.013 and 0.125). £ indicates the 

observations. 


TABLE 1 

Photographic Data for 29 Psc 


Plate 

Date 

Time (U.T.) 

Length (Min) 

D2061a 

Jan. 25, 1970 

01 h 46 m 

17 

D2062a 

Jan. 25, 1970 

02 h 15 m 

16 

D2062b 

Jan 25, 1970 

02 h 44 m 

24* 

D2074a 

Jan. 26, 1970 

01 h 47 m 

36* 


’’‘longer exposures due to variable clouds 
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TABLE 2 

Lines Found in 29 Psc 






*OBS 

X RMT 

Element 

; # 

\ 

X 0BS 

X RMT 

Element 

; # 

\ 

3691.52 

3691.557 

H 

18 


0.171 

4026.40 

4026.189 

He 

I 

(55)' 

10 . 673 

3696.88 

3697.154 

H 

17 


0.346 


4026.362 

He 

I 

(55). 

f 

3703.76 

3703.855 

H 

16 


0.617 

4076.4 

4076.0 

C 

II 

(36) 

0.034 

3711.92 

3711.973 

H 

15 


0.980 

4101.76 

4101.737 

H 

06 


7.937 

3721.72 

3721.940 

H 

14 


1.736 

4121.90 

4120.821 

He 

I 

(16)' 

10.089 

3734.14 

3734.370 

H 

13 


2.512 


4120.993 

He 

I 

(16). 

r 

3736.80 

3736.901 

Ca 

II 

( 3) 


4122.72 

4122.638 

Fe 

II 

(28) 


3750.08 

3750.154 

H 

12 


3.477 

4128.08 

4128.053 

Si 

II 

( 3) 

0.135 

3759.36 

3759.29 

Ti 

II 

(13) 

0.084 

4130.95 

4130.884 

Si 

II 

( 3) 

0.132 

3761.28 

3761.320 

Ti 

II 

(13) 

0.039 

4143.93 

4143.759 

He 

I 

(53) 

0.271 

3770.72 


H 

II 


4.738 

4169.00 

4168.971 

He 

I 

(52) 

0.039 

3797.94 


H 

10 


5.958 

4173.44 

4173.537 

Fe 

II 

(27) 

0.073 

3819.70 

3819.606 

He 

I 

(22) 

10 . 260 


4173.510 

N 

II 

(50) 



3819.761 

He 

I 

(22) 

I 

4178.32 

4178.39 

V 

II 

(25) 


3835.36 

3835.386 

H 

09 


7.116 

4179.05 

4178.855 

Fe 

II 

(28) 

0.095 







4206.00 

4205.92 

Ti 

II 

( 3) 




3843.64 

Cr 

I 

(87) 


4233.29 

4233.167 

Fe 

II 

(27) 

0.064 

3845.04 

3845.21 

S 

II 

(23) 


4267.15 

4267.02 

C 

II 

( 6)' 

10.067 

3848.98 

3848.96 

Cr 

I 

(69) 

0.018 


4267.27 

C 

II 

( 6). 

r 

3853.84 

3853.657 

Si 

II 

( 1) 

0.072 

4303.2 

4303.166 

Fe 

II 

(27) 

0.066 

3856.08 

3856.021 

Si 

II 

( 1) 

0.173 

4340.48 

4340.468 

H 

05 


8.310 

3862.72 

3862.592 

Si 

II 

( 1) 

0.145 

4351.53 

4351.764 

Fe 

II 

(27) 

0.032 

3867.68 

3867.477 

He 

I 

(20)' 

10.070 

4385.4 

4385.38 

Fe 

II 

(27) 

0.125 


3867.631 

He 

I 

(20). 

J 

4387.6 

4387.928 

He 

I 

(51) 

0.234 

3870.08 

3807.06 

A1 

II 

(74) 

0.035 

4395.28 

4395.031 

Ti 

II 

(19) 

0.061 







4416.84 

4416.817 

Fe 

II 

(27) 

0.063 

3889.12 

3889.051 

H 

08 


7.754 

4437.66 

4437.549 

He 

I 

(50) 

0.046 

3894.0 

3894.035 

Cr 

II 

(23) 


4469.96 

4469.92 

He 

I 

(15) 


3900.56 

3900.546 

Ti 

II 

(34) 

0.038 

4471.52 

4471.477 

He 

I 

(14)' 

10.346 


3900.680 

A1 

II 

( 1) 



4471.688 

He 

I 

(14) 

J 

3905.60 

3905.52 

Si 

I 

( 3) 


4475.12 






3906.72 

3906.95 

S 

II 

( 3) 


4481.26 

4481.128 

Mg 

II 

( 4)' 

10 . 263 

3908.88 

3908.76 

Cr 

I 

(23) 



4481.327 

Mg 

II 

( 4) 

J 

3913.44 

3913.46 

Ti 

II 

(34) 

0.065 

4488.2 

4489.185 

Fe 

II 

(37) 


3919.12 

3918.977 

C 

II 

( 4) 

0.037 

4508.18 

4508.283 

Fe 

II 

(38) 

0.045 

3920.64 

3920.677 

C 

II 

( 4) 

0.148 

4515.36 

4515.237 

Fe 

II 

(37) 


3924 .96 






4522. S 

4522.634 

Fe 

TI 

(38) 

0.041 

3933.670 

3933.664 

Ca 

II 

( 1) 

0.152 

4549.46 

4549.467 

Fe 

I T 

(38) 

0.108 

3968.50 

3968.47 

Ca 

II 

( 1) 


4555.92 

4555.890 

Fe 

II 

(37) 

0.035 

3970.02 

3970.074 

H 

07 


8.016 

4583.76 

4583.829 

Fe 

II 

( '8) 

0.087 

3978.90 

3978.87 

C 

II 

(37) 


4629.36 

4629.336 

Fe 

II 

(37) 

0.082 

4002.84 

4002.940 

V 

II 

( 9) 


4713.5 

4713.143 

He 

I 

(12) 

0.041 

4009.6 

4009.27 

He 

I 

(55) 

0.249 


4713.71 

He 

I 

(12 







TABLE 
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N(He)/N(H) from W^Hel) 


X 

N(He)/N(H) 

4026 

0.036 

4121 

0.037 

4388 

0.028 

4471 

0.015 

4713 

0.010 


<N(He)/N(H)> - 0.027 ± 0.012 


A 

CO 

3 ? 

CM 

© 

O 

m 

<0 

pH 

CO 

V to 

• 

• 

• 

• 

• 

« 

CM 

o 

V 

o 

tp 

s 

CO ^ 

o 
1— 1 

w 05 

CO 

Tp 

o 

© © © © © 

CO C» CO fl 
CM O CM O CM 


H 

• • 
CM r-l 

• 

o 

• • • • • 
CM 05 in H CM 

9 • • • • 

© o © o o 

• ••••••••••• 

CMHCMOt*OOCMOM*OCO® 








N 


(0 


(0 

0 ) 

c 


J 

I 

TP 

w 


a 


e 5 


CrJ 


■P 

o 

r-( Sh 

a v 

•H J 3 

■p e 

|H 3 
9 


fl 
B O 
O M 
■P 
< 



/*S 

to to 

CO 

© CO CM CM © 

Tp © i-l © © 


CW 

CO © 

© 

N CO CO 

© © © © © 


w 

© o 

CM 

fl rH © rH fl 

O O O O O 


fC 

• • 

• 



m 


o o 

o 

O O O O O 

O O O O O 

^H 






(0 







» 



© 



1 



o o oo 

© 

Gj 

o 

cm »r 


© © CM CM CM 

© © TP © to 

c 

H 

© tp 

CO 

CM N O M 1 CO 

© © © TP tP 

< 

X 

• • 

• 




>- 

CM CM 

00 

O O O fl fl 

r-l fl CO O 

0 






o 






c 




© 


C 3 


© 


© ip © 05 O 


-o 


OO tp 

m 

CM CM O 05 ^ 

O O © © to 

c 

00 

CM CO 

m 

NHOMO 

CM C» © © © 

3 


• • 

• 



.Q 


o © 

© 

O O O CM CM 

NHOOO 

< 






Vi 

• rs 

t** rH 


© tji M< tT CM 

to tp TP O © 

O 

a > 

CM TP 

CO 

© © © © © 

O P» © CO fl 

<H 

• 0 

CO o 

© 

© © © © © 

© © O H fl 


w ^ 

• • 

• 



*D 


© 00 

oo 

© © © © © 

OCHHH 

0 


rH rH 





coNto^niONiCMs m 
(onutot'CioonoooH^ 
OOOOOOOOOr-lTpo 


©ooooooooooo 


^t^o^pt^t^inocDCDoooo 

rlHNNNNOOOOOO 


HTTONOOOOOTf Nf>CHO 

®©o©©o©©©©o© 

ONHOOHNNlOirnH 

oooooooooooo 


oooo ooooooo 


co©©©©©tooocooo®® 


NNNNNNNNNNNN 


c n u* 

© CM 

tn <h 

TP Tp 


°fc 

CM 

Tp 


°a 

CM 

Tp 


%4 

CM 


CM 

N 


O O 
a O 

CM CM 
N N 


°P 

Tp 

N 


O O 
CM CM 

N N 


% 


I 


I 


I 


I 


I I 


I 


I I 


N 

I 


°a Q 



Cm a o 

a a cm cm 

© © 

Q 

Q Q 

© © © 

Tp TP TP TJ. 

ata 

© 

© © 



© © 

© 

© © 

ca a a 

•Q £> £2 JD 


to © 


© 

torp©toO©©0©t>Tp© 

c- to 

N © 

O © © © to 

© © © Tp © 

fl©©©©©©©©©©© 

© © 

© H 

© o © © o 

N © © © TP 

»t*HH^00O0t)00a'B© 

© i> 

rH 

© © © o © 

© fl © O CO 

©fl©©«»©©©©©« 

i-i © 

© 

© © © © © 

©©©OH 

fl©0©toto©©«Tp©0 

© © 

Tp 

« © © fl fl 

M 1 * CO (J) (J) 

Tp©©©rHrH©©©©©© 

© Tp 

TP 

© © « TP TP 

© © TP © © 



© 


CO 


CO 05 TP 

r-l H CO 


to 

CM 


oo c*- 

CM CO 


00 

CO 


00 

a 


M 

CQ 


0 

CM 






_ 


} 








TABLE 5 


Number Fractions £ ^ and log N^ 


Element 


29 Psc 

B stars 
log N t 

Cosmic 

H 

0.975 

12.00 

12.00 

12.00 

He 

2.45 x 10" 2 

10.40 

11.21 

11.16 

C 

2.04 x 10" 3 

8.32 

8.30 

8.48 

Mg 

2.08 x 10 -6 

6.33 

7.95 

7.46 

Si 

1.77 x 10" 4 

8.26 

7.45 

7.47 

Ti 

3.70 x 10" 6 

6.58 


4.82 



| 





Fe 


8.50 x 10 
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7.94 


6.90 
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